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The average lifetime {t[H)) it takes for a randomly started trajectory to 
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r) ' land in a small region {H) on a chaotic attractor is studied. t{H) is an 

important issue for controlling chaos. We point out that if the region H is 
visited by a short periodic orbit, the lifetime t{H) strongly deviates from 
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expanding eigenvalue of the short periodic orbit visiting H. 
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Controlling chaos by stabilizing one of the many unstable periodic orbits embedded 
within a given chaotic attractor is attainable with small, time-dependent changes in an 
accessible system parameter [|1|-|^. The idea is to observe a typical trajectory of the uncon- 
trolled system for some transient time, until it falls sufficiently close to the desired periodic 
orbit, and then to activate the control mechanism. An important issue related to the utiliza- 
tion of this method is the average lifetime of chaotic transients that precede the controlled 
periodic motion [|l],||-^]. 

Suppose that the uncontrolled chaotic attractor A describes the asymptotic behavior of 
the dynamical system O : D ^ D, -DC R*", also referred to as the original system. Let 
C, = O'^(^) be a point on a particular unstable periodic orbit 0J^ that we wish to stabilize. 
Furthermore, let the vicinity of the orbit H = H^{C,) be an ?7i- dimensional ball of radius 
e << 1 centered at ^. The probability that a randomly started trajectory does not reach H 
up to time t is ~ e'^'^'^^^\ The average lifetime r = t(H) is strongly correlated with the 
visitation frequency of typical trajectories to the region H, which is described in terms of 
the naturally invariant measure (/iAr) contained within H - fi]\f{H) |I|J^,§|. Obviously, if a 
certain region on a given chaotic attractor is visited more frequently by typical trajectories, 
the average lifetime it takes for an orbit to land in that region will be smaller. In the present 
study we address the following question: What is the deviation of r from fiN^H)^^ as a 
function of C, and e? 

We will demonstrate the existence of bursts in the lifetimes, i.e. significant deviations of 
r from fi]\f{H)^^, which appear when the H region encompasses a point on a short periodic 
orbit. In contrast to the overall r ^ ^n{H)~^ behavior, at these exceptional positions, the 
lifetime r is considerably prolonged as compared to ^n{H)~^. As the length of the shortest 
cycle visiting H increases, the parameter of this deviation, S{H) = T/fij^{H)^^, decreases 
rapidly towards 1. We will introduce a formula that relates the parameter S{H) to the 
repelling properties (expanding eigenvalue) of the shortest cycle within H. Furthermore, we 
will demonstrate that S{H) is independent of e (for e << 1). This is consistent with the 
previously reported scaling r ~ /iAr(i7)"^(see e.g. |jl]J 



The present study is motivated by the previous investigation of the logistic map with a 
hole 0. In this paper we present a theoretical explanation for the phenomenological result 
reported in Ref. [^] and generalize it to ID noninvertible and 2D invertible chaotic maps. 
(From our considerations a conjecture follows that similar phenomena occur generally in 
chaotic systems.) 

It will be useful to define an auxiliary modified map P,p!o[ 



, C'(e'), C e D\H 

M{0 = { \ (1) 

outside of the basin of A, ^' G H. 

A typical trajectory of the map O remains on the chaotic attractor forever, while the same 
trajectory in the map M eventually escapes through the region H, from now on also referred 
to as the hole. The average lifetime of chaotic transients created by the map M is equal to 
t{H), which we have defined above. Similar maps with a forbidden gap region arise in the 



context of communicating with chaos |Tl[] , and in calculation of the topological entropy [jT2 . 

To illustrate the concept of bursts, we consider two chaotic ID maps: (i) the asymmetric 
tent map 0{x) = kix, x < ki; 0{x) = k2{l — x), x > ki, ki,k2 > 0, k^-^ + A;^"*^ = 1, and 
(ii) the sinusoidal map 0{x) = sinvrx. Fig. |I] displays r as a function of the position of the 
hole C, {H = (C — e, ^ + e)), for the two paradigmatic maps (see also Ref. P). The width of 
the hole is kept constant (e = 0.005). Both graphs exhibit some common features: (i) the 
overall behavior of lifetimes follows the ^n{H)^^ pattern; (ii) strong local deviations from 
the ^]\f{H)~^ behavior - the bursts, observed as leaps in the lifetimes, occur when the hole 
interval sweeps across a short periodic orbit; (iii) the bursts are more significant if the length 
of the short periodic orbit is smaller. 

The explanation of the burst phenomenon requires the comparison of two concepts: 



i) the conditionally invariant measure |T3|-[T5[] (also referred to as the c-measure) - the 



concept associated with the modified system, and (ii) the naturally invariant measure |]^,|^ 
of the original system. In order to define these measures, imagine that we cover the chaotic 
attractor with cells (J) from a very fine grid. Then we randomly distribute a large number 
(A^) of points on the grid, and evolve them under the dynamics O for a long time T. Suppose 



that all initial points are colored blue, and that a point irretrievably changes color from blue 
to red immediately after its first entrance into the region H. Thus, the point xt = 0^{xq) at 
time T is blue if 0^{xo) ^ H for t G {0,1,..., T— 1}, and red otherwise. In the limit T — > cxd, 
the fraction of points found in a given cell /, is just the natural measure contained within 
that cell: /i^(/) = A/'t(/)/A^; A^ = EiNril), where Nril) = hT{I)+rT{I) denotes the total 
number of points in a given cell / at time T [§]. The number of blue (red) points within a 
given cell / at time T is denoted as &t(-^) {''^t{I))- The points which change color from blue 
to red under the action of the map O, are those which would escape the attractor under 
the action of the map M. Hence, if we were to evolve exactly the same initial conditions 
using the map M for the same time T, the number of surviving points (blue points) would 
be Bt = J2i bri,!) ~ A^exp (—T/r). In the limit A^ — * oo, T — > oo, the distribution of blue 
points converges to the c-measure of the modified system [p!3|-[T5[| . The fraction of blue points 



in a given cell / is simply the c-measure contained within that cell: ^c{I) = bT{I)/BT- 

The blue point at any time t > 0, was certainly not in the hole H at time t — 1. Therefore, 
M-\I) = 0-\I)\H and 6t(/) = &t-i(M"1(/)), which divided by Bt = fir-iexp (-l/r) 
yields the well known relation for the c-measure [^^-|TB 



ficil) = e^fic{M-\l)). (2) 



By summing the equation above over all the cells / we obtain 
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\n[l -i2c{H)]c:,fic{H). (3) 



We emphasize the importance of this observation. The average lifetime it takes for a typical 
trajectory to reach the small region H on the attractor is an inverse of the c-measure 
contained within that region {^c{H)^^) ■, which may significantly differ from the inverse of 
the natural measure {fij\f{H)^^). 

As an illustration, in Fig. |^ we display the c-measure for the modified version of the map 
sinvrx, in comparison to the natural measure of the original map. For the c-measure in Fig. 
^ a), the hole has been positioned at an arbitrary point, but not on the short periodic orbit. 



In this case, we observe that fic{H) ~ fiN^H), i.e., hn^H)'^ is a good approximation for 
the hfetime. In contrast, in Fig. |^ b) we display /i^ for the modified map sinvra;, with the 
hole positioned on the fixed point. We notice that iic{.H) strongly deviates from fiN^H). 
This case corresponds to the burst labeled 1 in Fig. [1| b). The overall agreement of the two 
measures is evident in both Figures ^ a) and b). However, at locations above the first few 
images of the hole, fic takes the shape of a well, with values which are considerably lower 
than fiN. When the hole lies on the fixed point (Fig. |^ b)), the wells are just above the 
hole itself. This results in a pronounced deviation of r = fic{H)~^ from fj,]\f{H)~^, which 
manifests as a burst. 

Now we compare the two measures globally. A chaotic repeller is a set of points on the 



attractor that never visit the hole |[Tj,|T5|JTT| . A trajectory that starts close to the repeller, 
does not escape the attractor for a long time. Therefore, the blue points at a large time T are 
located along the unstable manifold of the repeller. Their distribution along this manifold 



defines the c-measure [|T^ . Thus, the natural measure is constructed from all the points (at 
time T) on all the unstable manifolds, whereas the c-measure results only from points on 
parts of these manifolds. The parts which extend from the repeller up to the hole. As we 
reduce the size of the hole e, the repeller and its unstable manifold grow. Consequently, 
^c gradually approaches /iat, and for sufficiently small e, the two measures are practically 
identical. (For e = 0, fic becomes fij\[ [|^-|15|.) 

However, the deviation of the lifetime r = ^c{H)^^ from fi^^H)^^ depends only on the 
values of the two measures within the hole, and therefore is a local quantity. In order to 
make a more accurate comparison of fic and ^n, we introduce the following definitions. 
Consider a set P C -D such that /iAf(P) > 0. We define the quantity 

a{P)=ficiP)/MP), (4) 

which describes the relation between fiQ and fi^ within P. We also define the infiuence i{P) 
of the hole on the set P as 
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1{P) denotes the smallest integer for which the section 0~^^^\P) fl H becomes nonempty. 
The natural invariant measure within 0~^^^\P) is mapped to P in 1{P) iterates. The 
influence is just a fraction (0 < i{P) < 1) of ^in{P) that is mapped from the hole in the last 
1{P) time steps. 

Let Pe = Pe{S) C -D be an ?7i-dimensional ball of radius e (the same radius as the hole) 
centered at x. We ask the following question: Given a chaotic attractor and choosing the 
hole region, what is the behavior of a{P^) = fic{Pe)/fJ'N{Pe) as the position of P^ on the 
attractor is changed? 

By using Eq. (^ and the identity fi^i^P^) = /iAr(0^^(Pe)), we can write 

a{P,) = e^ ■ (1 - z(P,)) ■ a{M~'{P,)), (6) 

where / = /(P^) (in what follows, / = ^(-Pe))- 

Concerning the first factor in Eq. (^, note that the average lifetime typically scales 
like r ~ 1/6-^^(5) (^Dp{^) denotes the pointwise dimension at ^) [|T],^^, whereas the minimal 
number of iterates / for which 0^\P^) (1 H ^ ^ scales like / ~ ln(l/e) ||^. Therefore, 
exp(//r) ^ 1 + 1/t ~ 1. 

If the influence i{Pe) is small, the second factor in Eq. (^ is ^ 1. We argue that 
i{Pe), the influence of the hole on the region P^ decreases exponentially with /. For the 2D 
original map, O'^P^) is a narrow region which is stretched along the stable direction and 
squeezed along the unstable one 0]. The intersection of O^^P^) with the hole H is roughly 
a rectangle of length e and width eexp(— Ai/). For the ID map, 0^\Pe) fl if is an interval 
of width ~ eexp(— Ai/). In both cases, Ai denotes the positive Lyapunov exponent obtained 
for typical initial conditions on the attractor. Since the natural measure is concentrated 
along the unstable manifolds [0,^, we can relate yUAr((9^'(Pe) fl if ) ~ exp(— Ai/). Thus, due 
to the chaoticity of the map O we obtain i(Pe) ~ exp(— Ai/). 

Concerning the third factor in Eq. (^), we consider the set M~\P^) and the value 
a{M~'\P^)) in dependence of/. For the 2D maps, the set M~\P^) = 0^''{P^)\H is stretched 
exponentially fast with increasing / along the stable manifolds, and thus crosses many of the 
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unstable manifolds that carry both the natural and the c-nieasure. For the ID maps, the 
number of disjoint intervals that make the / — th preimage of P^ grows exponentially with /. 
Furthermore, they are scattered all over the attractor. Due to the chaoticity of the map O, 
in both cases M~\P^) becomes more democratic with larger Z, in the sense that the value 
a{M~\P^)) reflects the global agreement between the two measures. Thus, insofar as / is 
not small, a{M-\P^)) ~ 1. 

We conclude that for / = l{Pt) larger than some critical value (call it /c), all of the three 
factors in Eq. (^ are ~ 1 and therefore iic{Pe) — fJ'N{Pe)- This is consistent with the global 
agreement between the two measures. The critical value Ic depends on the chaoticity of the 
original map. For example, we may take Ic to be the smallest integer for which e^^^''" < 0.1 
(e.g. for the sin(7rx) map this gives /c ~ 3 — 4). When H maps to P^ in just a few iterates, 
so that / < /c, a significant difference is observed between Hc^Pe) and hn^P^) (this explains 
the wells in Fig. |^). 

Coming back to the average lifetimes, if the hole does not map back to itself in just a 
few iterates, i.e., if the shortest periodic orbit within H has a period larger than Z^, then 
fic{H) ^ fii\f{H), or simply r ~ fii\i{H)~^. This explains the overall behavior of lifetimes (see 
Fig. [I|). On the other hand, if H encompasses a short periodic orbit (period= 1{H) < l^), 
the two measures differ within the hole. Quantitatively, we substitute P^ -^ H in Eq. (^ 
and approximate e'^^^/"^ ~ 1 and /ic(M~'(-'^^) ^ fj.j\f{M~'-^^'>). This results in 

r^(l-^(/f))-V^(i7)-i. (7) 

Although 1{H) is small, the approximation a(M~'(^^) ^ 1 is justified if 1{M~'-^^^) > l^, or 
simply, if the period of the second shortest orbit within H exceeds Ic- We have tested relation 
(0) and consequently the approximations that lead to it in a number of systems. We have 
compared t{H) with fi]\f{H)^^ by changing the position of H from "the most exceptional" 
point, the fixed point, to longer cycles. In Fig. ^ we display a test of Eq. (|^ for the the 
generalized baker's map (see Ref. 0, p. 75, A^ = 0.35, \b = 0.40, a = 0.40, /? = 0.60), and for 
the Henon map (see Ref. |]16| a = 1.4,6 = 0.3). Recalling that i{H) decreases exponentially 



with 1{H), and considering Eq. (|^, we see that the parameter S{H) = {l—i{H))~^ decreases 
rapidly towards 1 with the increase of 1{H) (see Fig. ^). Eq. (0) is robust and can be apphed 
for holes of different shapes, as long as r >> 1. If (for the 2D maps) we tailor the hole as a 
rectangle with sides of length e parallel to the stable and unstable manifold segments, and 
center it on a short periodic orbit, then 

r^(l-A-i)-V^(iJ)-^ (8) 

A„ denotes the magnitude of the expanding eigenvalue of that orbit. Eq. (||) also applies 
to ID maps. Note that the approximation i{H) c^ A^^ assumes that the natural measure is 
smooth along the unstable direction within H. We observe that T/fij\f{H)~^ is independent 
of e. This is in accordance with the statement that the lifetime r scales with e just like 

Let us consider an application of Eq. (||). Suppose that we wish to control a chaotic 
system around an unstable fixed point. In order to obtain the position of the fixed point, 
its unstable eigenvalue, and other information required for the control, an observation of 
the free running system is needed |I|,0]. From this observation, we can also evaluate the 
visitation frequency to the e-vicinity of the fixed point, i.e. f^N^H = H^{^)). e is determined 
by the maximally allowed deviation of the control parameter from its nominal value [|l],3- 
The question of interest is how many iterates are needed on the average (r), before a chaotic 
trajectory enters the region H, when the control becomes attainable [0. The prediction 
given by fi]\f{H)^^ is an underestimate, since we are on the fixed point. For example, if the 
underlying dynamics of the system is the asymmetric tent map (with the same parameters 
as in Fig. |l|), and ii H = (0.62111801 . . .-0.002, 0.62111801 . . . + 0.002), the estimate for the 
lifetime ^^{H)^^ gives 250 iterates. On the other hand, the numerically calculated lifetime 
is ^ 627 iterates, which is more than twice as long. The lifetime obtained from formula (^ 
is 641 iterates, which is very close to the numerically calculated lifetime. Thus, Eq. (^) can 
be utilized to easily and accurately obtain r from an observation of the free running system. 

In summary, we have studied the average lifetime (r) it takes for a randomly started 
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orbit to land in a small region (H) on a chaotic attractor. That problem was introduced in 
Ref. [|I| as an important issue for controlling chaos. Our main result is that if a low-period 
unstable periodic orbit visits the region H, then the lifetime r significantly deviates from 
the inverse of the natural measure contained within H (/iAr(/7)~^). The parameter of this 
deviation, T/fj,N{H)~^, is a function of the expanding eigenvalue of that low-period orbit. 



FIGURES 
FIG. 1. t(^) (solid line) and jii<[{H)^^ (dashed line) vs. ^ for a) the asymmetric tent map 
(fcf 1 = 0.39, k^^ = 0.61), and b) the sinvrx map. H = {e. - e,e. + e), e = 0.005. 

FIG. 2. Figure displays ^c (solid line) in comparison to ^n (dashed line) for the sinvrx map 
(e = 0.005). Fig. 2a) shows ^c for the hole positioned on ^ = 0.66. The hole is not visited by a 
short periodic orbit. Note that ndH) ~ iJ,]\f{H). Fig. 2b) shows fic for the hole located on the 
fixed point. Note that ^c{H) < fiN{H). In both figures ^c and ^n are globally identical, except 
at the first 3-4 images of the hole, which are plotted underneath the graphs. 

FIG. 3. Numerically evaluated parameter t/ ixn{H)^^ for the Henon map (diamonds) and the 
Baker map (circles) in comparison to (1 — i{H))^^ (horizontal bars). The hole of radius e = 0.005 
is centered on the shortest periodic orbit (period= 1(H)) visiting H. 
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